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• We propose a new technique, called quantum optical coherence tomogra- 

o 

■ phy (QOCT), for carrying out tomographic measurements with dispersion- 

cancelled resolution. The technique can also be used to extract the frequency- 
dependent refractive index of the medium. QOCT makes use of a two-photon 
^ • interferometer in which a swept delay permits a coincidence interferogram 

pH . to be traced. The technique bears a resemblance to classical optical coher- 

ed . 

^ . ence tomography (OCT). However, it makes use of a nonclassical entangled 

twin-photon light source that permits measurements to be made at depths 
greater than those accessible via OCT, which suffers from the deleterious ef- 
fects of sample dispersion. Aside from the dispersion cancellation, QOCT 
offers higher sensitivity than OCT as well as an enhancement of resolution by 
a factor of 2 for the same source bandwidth. QOCT and OCT are compared 
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using an idealized sample. 
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I. INTRODUCTION 



Optical Coherence Tomography (OCT) has become a versatile and useful tool in biophoton- 
ics IQ . It is a form of range- finding that makes use of the second-order coherence properties 
of a classical optical source to effectively section a reflective sample with a resolution 
governed by the coherence length of the source. OCT therefore makes use of sources of 
short coherence length (and consequently broad spectrum), such as superluminous LEDs 
and ultrashort-pulsed lasers. 

A number of non-classical (quantum) sources of light have been developed over the past 
several decades |^ and it is natural to inquire whether making use of any of these might be 
advantageous. The answer turns out to be in the affirmative. Spontaneous parametric down- 
conversion (SPDC) [Q is a nonlinear process that generates entangled beams of light; these 
have been utilized to demonstrate a number of two-photon interference effects that cannot 
be observed using traditional light sources 0. We demonstrate that such entangled-photon 
fourth-order interference effects may be used to carry out range measurements similar to 
those currently obtained using classical OCT, but with the added advantage of even-order 
dispersion cancellation 0. This is possible by virtue of the non-classical nature of the 
light produced by SPDC. We refer to this new technique as quantum optical coherence 
tomography (QOCT). 

II. CLASSICAL OPTICAL COHERENCE TOMOGRAPHY (OCT) 

The sample investigated in the course of our calculations, classical and quantum alike, is 
represented by a transfer function H . This quantity describes the overall reflection from all 
structures that comprise the sample. For an incident undepleted monochromatic plane wave 
of angular frequency cj, 

H{uj)= / dzr{z,uj)e'^^^''''\ (1) 
Jo 
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Here r{z,uj) is the complex reflection coefficient from depth z and (p{z,uj) is the phase 
accumulated by the wave while travelling through the sample to the depth z. 

The basic scheme of OCT [|1| is illustrated in Fig. 1. We assume that the classical 
source produces cw incoherent light with a short coherence time of the order of the inverse 
of its spectral width (the results also apply to the case of a pulsed source, however). We 
characterize the source S with a power spectral density S{ujq + Q) where ujq is its central 
angular frequency. The light is divided by a beam splitter into the two arms of a Michelson 
interferometer. A variable delay r, imparted by a scanning mirror, is placed in the "reference 
arm" while the sample is placed in the "sample arm" . The reflected beams are recombined 
by the beam splitter and an interferogram /(r) is measured: 

/(r) oc To + 2Re{r(r)e-''"''"}. (2) 

The self interference term Fq and the cross- interference term F(r) are given by 

Fo = ldn[i + \H{ujo + s{n), (3) 

and 

F(r) = JdQ H{uJo + Q) S{Q) e~'^^ = /i,(r) * s(r), (4) 

respectively, where hc{T) is the inverse Fourier transform of H{ujq + Q) with respect to fl, 
and s(r) is the correlation function of the source [the inverse Fourier transform of S'(n)]. 
The symbol * represents the convolution operation. 

The physical underpinnings of this scheme may be understood by examining the inter- 
ference of light propagating in the two paths created by the beam splitter (Fig. 1). A 
monochromatic wave of frequency ujq + Q emitted from S acquires a reflection coefficient 
H(ujq + fi) in the sample arm, but only a phase factor in the reference arm. As is 

clear from Eq. (Q), the resulting interferogram includes a self- interference contribution from 
the two paths (Eq. ^): a factor of unity from the reference path and a factor of \H{ujq + P 
from the sample path. The cross- interference contribution, which resides in Eq. (^, is 
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the product of these two terms, H{ujq + Q) and e'*^'^''"'"^)'^ (one is conjugated, but this is of 
no significance in OCT). This term may also be expressed as a convolution of the sample 
reflection with the coherence function of the source, the temporal width of which serves to 
limit the resolution of OCT. 

III. QUANTUM OPTICAL COHERENCE TOMOGRAPHY (QOCT) 

The scheme we propose for QOCT is illustrated in Fig. 2. The twin-photon source is 
characterized by a frequency-entangled state given by 

1^) = Jdn c{n) \ujo + n,ujo-n), (5) 

where Q is the angular frequency deviation about the central angular frequency uq of the 
twin-photon wave packet, C(^) is the spectral probability amplitude, and the spectral dis- 
tribution of the wave packet S{Q) = is normalized such that / dQS{Q) = 1. For 
simplicity, we assume is a symmetric function and that both photons reside in a common 
single spatial and polarization mode. 

Interferometry is implemented by making use of a seminal two-photon interference exper- 
iment, that of Hong, Ou, and Mandel (HOM) 0. The HOM beam-splitter interferometer is 
modified by placing a reflective sample in one of the paths in the interferometer and a tem- 
poral delay r is inserted in the other path, as shown in Fig. 2. The two photons, represented 
by beams 1 and 2, are then directed to the two input ports of a symmetric beam splitter. 
Beams 3 and 4, the outputs of the beam splitter, are directed to two single-photon-counting 
detectors, Di and D2. The coincidences of photons arriving at the two detectors are recorded 
within a time window determined by a coincidence circuit. The delay r is swept and the co- 
incidence rate C(r) is monitored. If a mirror were to replace the sample, sweeping the delay 
would trace out a dip in the coincidence rate whose minimum would occur at equal overall 
path lengths, which we define as zero delay. This dip is a result of quantum interference of 
the two photons within a pair. 
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For a sample described by H{u), as provided in Eq. (|l|), the coincidence rate C(r) is 
given by 

C(r) oc Ao-Re{A(2r)}, (6) 

where the self- interference term Aq and the cross- interference term A(r) are defined as 
follows: 

Ao = J dn\H{uo + Q)\^ S{Q), (7) 

and 

A(r) = JdQ H{uJo + Q) H*{uJo - Q) S{Q) e''^^ = hgir) * s(r). (8) 
Here hg{T) is the inverse Fourier transform of Hg{Q) = H{ujq + ^7) H*{ujq — Vt) with respect 

to n. 

It is important to highlight the distinctions and similarities between Eqs. (|^), and 
(|), and Eqs. (|^), (|^), and (^. The unity OCT background level in Eq. is, fortuitously, 
absent in Eq. (|^) for QOCT. Moreover, the QOCT cross-interference term in Eq. (^) is 
related to the reflection from the sample quadratically; the sample reflection is therefore 
simultaneously probed at two frequencies, ujq + VL and ujq — fi, in a multiplicative fashion. 
Finally, the factor of 2 by which the delay in the QOCT cross-interference term in Eq. 
is scaled, in comparison to that in Eq. (0) for OCT, leads to an enhancement of resolution 
in the former. 

A particularly convenient twin-photon source makes use of spontaneous parametric down- 
conversion (SPDC) [^. This process operates as follows: a monochromatic laser beam of 
angular frequency Up, serving as the pump, is sent to a second-order nonlinear optical crystal 
(NLC). Some of the pump photons disintegrate into pairs of downconverted photons. We 
direct our attention to the case in which the photons of the pairs are emitted in selected 
different directions (the non-collinear configuration). Although each of the emitted photons 
in its own right has a broad spectrum, by virtue of energy conservation the sum of the 
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frequencies must always equal uip. Because of the narrow spectral width of the sum frequency 
(which is the same as the pump frequency), the photons interfere in pairs. But because of 
the broadband nature of each of the photons individually, they serve as a distance-sensitive 
probe not unlike the broadband photons in conventional OCT. 

IV. COMPARISON OF QOCT AND OCT 

The sample model presented in Eq. (|l]) may be idealized by representing it as a discrete 
summation 

Hico) = Y.r,{co)e''^^^-\ (9) 
j 

where the summation index extends over the layers that constitute the sample. This is a 
suitable approximation for many biological samples that are naturally layered, as well as for 
other samples that are artificially layered such as semiconductor devices. This approximation 
is not essential to the development presented in this paper, however. 

We further assume, without loss of generality, that the dispersion profile of the media 
between all surfaces are identical, so that (Pj{uj) = (3{uj)zj, where = n{uj)uj/c is the 
wave number at angular frequency cj, Zj is the depth of the j*^ layer from the sample surface, 
n{uj) is the frequency-dependent refractive index, and c is the speed of light in vacuum. We 
expand /?(co'o + ^) to second order in Vt: (3{uJo + Q) ^ f3o + f3'Q + where is the 

inverse of the group velocity vq at uq, and P" represents group velocity dispersion (GVD) 
i- 

In the case of OCT, using Eqs. @ and leads to a cross-interference term given by 

^(r) = E-.4°t--2^)e"^"^ (10) 

where s^^''\-) arises from reflection from the j*'' layer after suffering sample GVD over a 
distance 2zj, the subscript d indicates dispersion, and the superscript (Oj) indicates that 
dispersion is included from the surface of the sample (0) all the way to the j^^ layer. The 
quantity s^"''^^(-) is thus the Fresnel transformation of S{Q) with dispersion coefficient (3" 0: 
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The effectiveness of OCT is tfierefore limited to samples tliat do not exliibit appreciable 
GVD over the depth of interest. 

In the case of QOCT, on the other hand, Eqs. (||) and (|) result in a cross-interference 
term given by the sum of two contributions: 



the first contribution represents reflections from each layer without GVD, while the second 
contribution represents cross-terms arising from interference between reflections from each 
pair of layers. The quantity s(-) is the correlation function of the source defined previously, 
and the quantity s^J''\-) is the Fresnel transformation given in Eq. (pH]). In contrast to OCT, 
only dispersion between the j^^ and k^^ layers survives, as is evident by the superscript (jk). 
The terms comprising the first contribution in Eq. (|12D include the information that is often 
sought in OCT: characterization of the depth and reflectance of the layers that constitute 
the sample. The terms comprising the second contribution in Eq. (|1^) are dispersed due to 
propagation through the inter-layer distances Zj—Zk] however, they carry further information 
about the sample that is inaccessible via OCT. Two complementary approaches can be used 



to extract information from Eq. (^21) : 1) averaging the terms that comprise the second 
contribution by varying the pump frequency while registering photon coincidences such that 
the exponential function averages to zero, which leads to unambiguous optical sectioning 
information resident in the first contribution; and 2) isolating and identifying the terms of 
the second contribution to obtain a more detailed description of the sample than is possible 
with OCT. This can, in fact, be achieved by making use of the Wigner distribution as will 
be demonstrated at the end of this section. 

We now proceed to provide a numerical comparison between QOCT and OCT using Eqs. 
(p!0|) and (|l^). Consider a sample comprising of two reflective layers buried at some depth 
below the surface of a medium, as illustrated at the very top of Fig. 3. For the purposes of 
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our calculation, we arbitrarily choose amplitude reflection coefficients ri = 0.1 and r2 = 0.2, 
separation distance di = 10 /im, and depth below the sample surface = 0.1 mm. For both 
OCT and QOCT, calculations are carried out by assuming that the source has a central 
wavelength Aq = 2ttc/u!o = 812 nm and a Gaussian spectral distribution with a bandwidth 
(full width at 1/e of maximum) of 155 nm, which corresponds to a wave packet of temporal 
width 14 fsec and length 4.2 fim in free space. In the context of QOCT, this can be realized 
by means of a /5-barium borate NLC of thickness 1 mm pumped by a source of wavelength 
Xp = 406 nm. Using type-I SPDC, a NLC cut at an angle 29° with respect to the optic axis 
generates light centered about the degenerate wavelength Aq = 2Ap = 812 nm. For purposes 
of illustration, we neglect reflection from the top surface of the sample and assume that the 
sample dispersion profile is characterized by: (3' = 5x 10^^ sm^^ and (3" = 1.8 x 10^^^ s^ m^^. 
These correspond to a highly dispersive material; however, common materials such as "higher 
dispersion crown glass" are an order-of-magnitude more dispersive. 



The results of this calculation are displayed in Fig. 3 for OCT (thin rapidly varying 
gray curve) and QOCT (black broken curve representing the full signal; black solid curve 
representing the signal averaged over pump frequency). Because of dispersion it is clear 
that no useful information about the sample is available from OCT. QOCT, on the other 
hand, yields a pair of high-resolution dispersion-cancelled coincidence-rate dips at delays 
corresponding to reflections from the two surfaces. Moreover, the QOCT resolution is a 
factor of 2 superior to that achievable via OCT in a dispersionless medium. The peak 
between the two dips evidenced in the full QOCT signal (black broken curve), which could 
alternatively be a dip depending on the phases of the terms in the second contribution in Eq. 
(|l^), is a result of quantum interference between the probability amplitudes arising from 
reflection from the two different surfaces. This is in contrast to the black solid-curve dips, 
which are a result of quantum interference between the probability amplitudes arising from 
reflection from each surface independently. The breadth of the middle peak is determined 
only by the dispersion of the medium residing between the two reflective surfaces and not 
by the nature of the material under which they are buried. It is clear, therefore, that the 
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dispersion of the region between the two surfaces may be determined by measuring the 
broadening of the middle peak in comparison with the two dips. 

It is worthy of note that dispersion cancellation occurs for all even powers of the expansion 
of l3{uj). Thus if the phases of reflection from the surfaces are random, which provides a 
model for transmission through a turbid or turbulent medium, only the middle peak will 
wash out, while the dips arising from reflections from the surfaces of interest are unaffected. 
In OCT, such random phase variations serve to deteriorate, and possibly destroy, information 
about the sample. 

In Fig. 4 we plot results for the same example examined above, except that one of the 
layers of interest is situated at the surface of the sample rather than being buried beneath 
it. In this case OCT gives intelligible results although the return from the second layer is 
clearly broadened as a result of dispersion. On the other hand, the results for QOCT are 
identical to those shown in Fig. 3 for the same two-layer object buried under a dispersive 
medium. QOCT is also seen to exhibit higher sensitivity than OCT for weakly reflective 
samples. This is because the self-interference term in QOCT [Eq. (|^] does not include the 
factor of unity present in the self-interference term of OCT [Eq. (|^)]. 

Finally, we address the use of the Wigner distribution for extracting information about 
the sample via the QOCT cross- interference term A. Examining Eq. (|^), and assuming that 
the bandwidth of S{Q) is greater than that of Hg{Q), we obtain 



This is precisely the Wigner distribution function of the function H{Q), with parameters 



1 of the interferometer in Fig. 2. The quantity ujq may be changed by varying the pump 
frequency Up = 2uJo. Although this technique might be expected to face practical difficulties 
because the direction of SPDC changes as the frequency is varied, this could be mitigated 
by adjoining a wave- guiding mechanism to the twin-photon source, as is customary when 




(13) 



T and ujq |jTT|. Knowledge of A(r, cuq) for all relevant values of r and uq guarantees that 



H{Q) may be reconstructed |T^. The quantity r is varied by changing the delay in path 
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using periodically-poled NLCs, for example. Furthermore, such an approach would enable 
the output light to be directly coupled into an optical fiber and thus integrated into systems 
already familiar to the practitioners of OCT. 



V. ADVANCED- WAVE INTERPRETATION 

The operation of QOCT may be understood in a heuristic way by considering an 
advanced-wave interpretation similar to that employed by Klyshko in the context of spatial 
interferometers [|l^]. In such an interpretation one of the detectors may be thought of as 
being replaced by a classical light source with its waves traced backward through the optical 
system, and the twin-photon source may be thought of as a refiector. The intensity measured 



at the location of the other detector then mimics the coincidence rate [jT3|. Applying this 
interpretation to QOCT, assume that Di in Fig. 2 is replaced by a classical light source that 
emits a monochromatic wave of frequency u^o + ^- The beam splitter results in this wave 
being partitioned into two paths. In one of these (path 1) the wave travels backward through 
the delay, changes direction and fiips its frequency about ujq to cuq — ^ at the twin-photon 
source, and then propagates forward through path 2. Finally, it refiects from the sample 
and reaches D2 having acquired a weig hting factor of e*('^o+^)^ H{uJo - fi). The second wave 
(path 2 after the beam splitter) reflects from the sample, changes direction and frequency 
from tuo + n to cjq ~ ^ at the source, and then undergoes a delay r in path 1 en route to D2, 
acquiring a weighting factor e^^^°~^^'^ H{ujq + Vt). The self-interference contribution in Eq. 
(0) is given by the sum of the squared amplitudes of these two terms. The cross-interference 
contribution to the interferogram is one of these terms multiplied by the complex conjugate 
of the other: e"*^^"^ H{uq + Vt) H*{ijJq — Vt). This interpretation makes clear the origin of the 
salutary time-scaling by a factor of 2 and the absence of interference fringes at frequency ujq 
from the QOCT interferogram. Both QOCT interfering waves reflect from the sample and 
they do so at conjugate frequencies, whereas in OCT one of the waves reflects from a mirror, 
which gives rise to the deleterious unity term that is absent from QOCT. In contrast to 
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OCT, the complex conjugate present in the cross-interference term is of central importance 
in QOCT. 



VI. CONCLUSION 

We have presented a new technique, called quantum optical coherence tomography (QOCT), 
which utilizes the wave-packet nature of photons generated in pairs via spontaneous paramet- 
ric down-conversion (SPDC). Each photon of the pair inherently occupies a broad spectrum 
even though the pump is monochromatic: the bandwidth is determined by the length of 
the nonlinear crystal. QOCT yields performance superior to that of a classical optical co- 
herence tomography (OCT) on three counts: 1) the resolution is enhanced by a factor of 2 
for the same source bandwidth; 2) it has greater sensitivity for weakly reflecting samples; 
and 3) sample group-velocity dispersion does not result in a deterioration of resolution with 
increasing depth into the sample. Moreover, the frequency dependent refractive index of the 
medium, which is inaccessible to OCT, may be extracted. 
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and Imaging Systems (CenSSIS), an NSF engineering research center. 
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FIGURE 1 



FIG. 1. Setup for optical coherence tomography (OCT). BS stands for beam sphtter; D is a 
detector; and r is a temporal delay introduced by moving the reference mirror. 
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FIG. 2. Setup for quantum optical coherence tomography (QOCT). BS stands for beam splitter 
and r is a temporal delay. Di and D2 are single-photon-counting detectors that feed a coincidence 
circuit. 
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FIG. 3. Normalized intensity /(r) (thin rapidly varying gray curve; left ordinate) and normal- 
ized coincidence rate C(r) (thick black curves; right ordinate) versus normalized delay (scaled by 
half the group velocity vq/2) for a two-layer sample buried under a dispersive medium. The black 



broken curve represents the full QOCT signal [Eq. (12)] whereas the black solid curve represents 
the QOCT signal after averaging over the pump frequency [Eq. (|l2[), first contribution]. The black 
broken curve coincides with the black solid curve everywhere except where the black broken curve 
is visible. The structure of the sample is shown at the top of the figure. The OCT signal yields no 
useful information, whereas the QOCT signal, by virtue of the dispersion-cancellation properties 
of this technique, clearly reveals the presence of the surfaces in the sample. 
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FIG. 4. Normalized intensity /(r) (left ordinate) and normalized coincidence rate C(t) (right 
ordinate) versus normalized delay for a two layer sample at the surface of a medium. Curves have 
the same significance as in Fig. 3. 
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